Superconductivity in the doped topological insulator Cu2;Bi2Se3 under high pressure 
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We report a high-pressure single crystal study of the topological superconductor Cui,Bi2Se3. 
Resistivity measurements under pressure show superconductivity is depressed smoothly. At the same 
time the metallic behavior is gradually lost. The upper critical field data BciiT) under pressure 
collapse onto a universal curve. The absence of Pauli limiting and the comparison of Bc2{T) to 
a polar state function point to spin-triplet superconductivity, but an anisotropic spin-singlet state 
cannot be discarded completely. 
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Recently, topological insulators have sparked a wide 
research interest, because they offer a new playground for 
the realization of novel states of quantum matter P, Q . 
In 3D topological insulators (TI) the bulk is insulating, 
but the 2D surface states — protected by a non-trivial Z2 
topology — are conducting. Exemplary TIs are the well- 
known thermoelectric effect materials Bi2Te3 and Bi2Se3, 
where photoemission [3, 4] and magnetotransport Q 
experiments have demonstrated the presence of a Dirac 
cone in the energy dispersion of the surface states, as in 
graphene, and 2D signatures like the quantum Hall effect. 

The concept of topological insulators can also be ap- 
plied to superconductors (SCs), due to the direct analogy 
between topological band theory and superconductivity: 
the Bogoliubov - de Gennes Hamiltonian for the quasi- 
particles of a SC has a close similarity to the Hamilto- 
nian of a band insulator, where the SC gap corresponds 
to the gap of the band insulator 0,11. All topological 
SCs in 1, 2 and 3D are predicted to be non-trivial SCs 
with odd-parity Cooper pair states 0, [l3|. Of major 
interest in the field of topological SCs is the realization 
of Majorana fermions P, that are predicted to exist 
as protected bound states on the edge of the 1, 2 or 3D 
superconductor. Majorana fermions are of great poten- 
tial interest for topological quantum computation P, Q ■ 
However, topological SCs are scarce. The B phase of 
•^He has recently been identified as an odd-parity time- 
reversal invariant topological superfluid (class 3D Z) , 
whereas the correlated metal Sr2Ru04 is a time-reversal 
symmetry breaking chiral 2D p-wave SC (class D Z) Q. 
Here we focus on a new potential candidate, namely the 
doped topological insulator Cui:Bi2Se3 fl^ . 

Recently, Hor and co-workers initiated a new route to 
fabricate topological superconductors, namely by react- 



low Tc, which cast some doubt on the bulk nature of SC. 
As regards Cua;Bi2Se3, this concern was taken away by 
Kriener ei a/. 14,15|. Bi2 Sea single crystals electrochem- 



ing the TIs Bi2Se3 and Bi2Te3 with Cu or Pd [12,ll3|. By 
intercalating Cu^"*" into the Van der Waals gaps between 
the Bi2Se3 layers, SC occurs with a transition temper- 
ature Tc = 3.8 K in Cu:,Bi2Se3 with 0.12 < x < 0.15. 
However, the reported SC shielding fractions were rather 
small and the resistance never attained a zero value be- 



ically intercalated by Cu^+ have a SC volume fraction of 
about 60 % (for x — 0.29) as evidenced by the large jump 
in the electronic specific heat at Tc- SC was found to be 
robust and present for 0.1 < x < 0.6. 

Photoemission experiments conducted to study the 
bulk and surface electron dynamics reveal that the 
topological character is preserved in Cu3;Bi2Se3 flif . 
Based on the topological invariants of the Fermi sur- 
face, Cu2:Bi2Se3 is expected to be a time-reversal invari- 
ant fully-gapped odd-parity topological SC 0, [13] ■ A 
recent study of the Cooper pairing symmetry within a 
two-orbital model led to the proposal that such a state 
can be favored by strong-spin orbit coupling (17[ . Indeed 
several experiments have revealed properties in line with 
topological SC. The magnetization in the SC state has 
an unusual field variation and shows curious relaxation 
phenomena, pointing to a spin-triplet vortex phase [3l ■ 
Point contact measurements reveal a zero-bias conduc- 
tance peak in the spectra, which possibly provides a first 
evidence of Majorana fermions [19[ . These signatures of 
topological SC make the experimental determination of 
the basic SC behavior of Cua;Bi2Se3 highly relevant. 

In this Letter we report the response of the supercon- 
ducting phase of Cuo.3Bi2.1Se3 to high pressure. SC is 
depressed smoothly. The T variation of the upper critical 
field Bc2{T) shows a universal behavior as a function of 
pressure, which lies above the standard variation for an 
s-wave SC. In addition, the absence of Pauli limiting and 
the Bc2{T) variation in agreement with the polar state 
model both point to spin-triplet SC. 

Single crystals of Cu2;Bi2Se3 were prepared by melt- 
ing high purity elements at 850 °C in sealed evacuated 
quartz tubes, followed by slowly cooling till 500-600 °C. 
After growth the samples were annealed for 60-100 hours. 
The best samples — with a zero resistance state below Tc 
— were obtained for a nominal Cu content x ~ 0.3, a 
nominal Bi content 2.1, and rapidly quenching after an- 
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FIG. 1; (color online) Resistivity of Cuo.3Bi2.1Se3 as a func- 
tion of temperature around Tc at pressures up to 2.31 GPa as 
indicated. Upper inset: p{T) at p = for sample #1. Lower 
inset, left axis: i?jv(293 K) (closed symbols) and Rn{^ K) 
(open symbols) as a function of pressure for samples #1 and 
#2. The subscript A'' means the data are normalized at 
i?(293 K, p = 0). Lower inset, right axis: (triangles) Car- 
rier concentration n(p) at T = 4 K for sample #3. 



nealing. A typical resistivity trace for a current in the ab- 
plane of the rhombohedral crystal structure {R3m space- 
group) is presented in the upper inset of Fig. 1. The re- 
sistivity shows metallic behavior, the residual resistance 
Pq — 0.15 mrjcm is relatively low, and the carrier concen- 
tration n ~ 1.2 X lO^*' cm^^. These transport parameters 
are similar to those reported in Refs. 
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The high-pressure measurements were carried out us- 
ing a hybrid clamp cell made of NiCrAl and CuBe al- 
loys. Samples were mounted on a plug which was placed 
in a Teflon cylinder with Daphne oil 7373 as hydro- 
static pressure transmitting medium. The pressure cell 
was attached to the cold plate of a ^He refrigerator 
{Tbase — 0.24 K). The ac-resistivity (/ — 13 Hz) was mea- 
sured, using the lock-in technique with a low excitation 
current (/ = 100 /iA), on two plate-like single-crystals 
with the current in the ab plane. The suppression of SC 
in a magnetic field was investigated by resistance mea- 
surements in fixed fields B \\ ab (sample #1) and B \\ c 
(sample ^2). The pressure was determined by measuring 
the SC transition temperature Tdp) of Pb in a separate 
experiment [20]. 

In Fig. 1 we show p{T) around Tc for sample ^1 at 
pressures up to 2.31 GPa. The normal state resistance 
i?(4 K) shows a pronounced increase as a function of 
pressure, while Tc steadily decreases. The overall good- 
quality of sample #1 is attested by the relatively small 
width ATc (as measured between 10 % and 90 % of 
i?(4 K)), which ranges from 0.25 K at p = to 0.06 K 
a.t p = 2.31 GPa. Nevertheless, a tail towards low tem- 
peratures associated with ^ 10 % of the resistance path 




FIG. 2: (color online) Superconducting transition tempera- 
ture as a function of pressure for samples #1 and ^2 as indi- 
cated. The solid line represents a polynomial fit with linear 
and quadratic terms which serves to extrapolate the data. 



is present. For sample #2 the resistance does not reach 
R = 0, but remains finite at the level of 10 % of i?(4 K). 
In Fig. 2 we report the variation Tc{p), determined by 
the midpoints of the resistive transitions, which almost 
coincide for both samples. The solid line in Fig. 2 (see 
caption) suggests Tc might be suppressed at a critical 
pressure pc as high as ~ 6.3 GPa. 

The depression of Tc can be understood qualitatively 
in a simple model for a low carrier density SC where 
Tc ez3exp[-l/A^(0)Vb], with Qd the Debye tempera- 
ture, iV(0) ^ m*'n}/^ the density of states (with m* the 
effective mass) and Vo the effective interaction param- 
eter The increase of i?(4 K) under pressure by a 
factor > 5 indicates a decrease of the carrier concentra- 
tion n, which in turn leads to a reduction of iV(0) and 
Tc- The reduction of n is also apparent in R{T) which 
gradually loses its metallic character (see lower inset in 
Fig 1): for p > 2.02 GPa i?(4 K) exceeds i?(293 K) for 
sample #2. Hall effect data taken at T = 4 i^T on a third 
sample (#3) for i? 1 1 c confirm n decreases under pressure 
(see lower inset in Fig. 1). 

Next we turn to the pressure variation of the upper 
critical field Bc2{T). The depression of SC by a magnetic 
field was measured for B || ab (i?°| — sample #1) and B \\ 
c {B^2 ~ sample #2). Tc{B) extracted as the midpoints 
of the transitions is reported for each pressure in Fig. 3. 
We first discuss the results at p — 0. The upper criti- 
cal field shows a moderate anisotropy for a layered com- 
pound with B^^{T ^ 0) = 5.6 T and B^^i^ 0) = 1.9 
T. The anisotropy parameter 7°" = B^^/B^^ = 2.9. Us- 
ing the relations B^^ = '^o/2n^lb and B^^ = $o/27r^a6^c, 
where $0 is the fiux quantum, we derive SC coherence 
lengths ^ab ~ 13 nm and = 4 nm. These values are 



in good agreement with those reported recently [1^ [14 1 . 
It is important to distinguish whether our samples are in 
the clean or dirty limit, or in between. Notably, because 
a sufficiently clean sample with an electron mean free 
path £ larger than ^ is a prerequisite for triplet Cooper 
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FIG. 3: (color online) Temperature variation of the upper 
critical field Bc2{T) for B \\ ab and B |j c at pressures of 0, 
0.26, 0.67, 1.42, 2.02 and 2.31 GPa (from top to bottom). 
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FIG. 4: (color online) Upper critical field Bc2{T) divided by 
Tc and normalized by the initial slope \dBc2/dT\T^ as a func- 
tion of the reduced temperature T/Tc at pressures of 0, 0.26, 
0.67, 1.42, 2.02 and 2.31 GPa for B \\ ab (closed symbols) and 
-B II c (open symbols). The red and blue full lines represent 
model calculations for an s- and p-wave superconductor. The 
black curve matches the data closely and represents the po- 
lar state model function scaled by a factor 0.95, which would 
result from a 5 % larger initial slope in the model. 



pairing [22] 
the relation 



An estimate for £ can be obtained from 
= hkp / pone'^ , assuming a spherical Fermi 
surface Sp = 4:TTkp with wave number kp — (STr^n)^/'^. 
With n = 1.2 X 10^^ m^^ and po = 1.5 x 10"^ nm (sec 
Fig. 1) we calculate kp — 1.5 x 10^ m~^ and £ — 3A nm, 
which ensures £ > S,. 

A more detailed analysis can be made employ- 
ing the slope of the upper critical field dBcz/dT at 
Tc (23j . In the dirty limit the initial slope is given by 
\dBc2/dT\T^ = 4480 -7^07 where 7 is the Sommerfeld co- 
efficient of the specific heat. With 7 = 22.9 J/m^K^ |1J] 
we calculate \dBc2/dT\T^ = 0.15 T/K. This value is 
much lower than the measured values 2.0 T/K {B \\ ah) 
and 0.6 T/K [B \\ c), which confirms our samples are 



not in the dirty limit. By adding the clean limit term 
\dBc2/dT\T, = 1.38 X 10^^ • -i'^TjSl in the model ^ 
estimates for £ and ^ can be extracted from the experi- 
mental values of \dBc2/ dT\T^. For B \\ ab (c) we obtain 
^ ~ 90 (45) nm and ^ ~ 9 (19) nm. Notice, in this analy- 
sis we used the normal-state 7-value as input parameter. 
If we consider that only part of the sample becomes su- 
perconducting, a reduced 7^ value [lj| should be used. 
This will affect the absolute values of the deduced pa- 
rameters, but not our conclusion £ > ^. We conclude our 
samples are sufficiently pure to allow for odd-parity SC. 

Under pressure Bc2{T) gradually decreases and the 
anisotropy parameter reduces from 7"" — 2.9 at p = to 
2.1 at the highest pressure p — 2.31 CPa. The coherence 
lengths increase to = 15 nm and S.c ~ 7 nm. As men- 
tioned above the increase of po and the gradual loss of 
metallic behavior under pressure can for the major part 
be attributed to a corresponding decrease of n, which 
tells us the ratio £/£_ > 1 in the entire pressure range. 

The functional behavior of Bc2 does not change with 
pressure, as is demonstrated in Fig. 4. All the Bc2{T) 
curves collapse on one single universal function b*{t), 
with b* = {Bc2/Tc)/\dBc2/dT\T, and t = T/T^ the re- 
duced temperature. This holds for B \\ ab, as well as 
for B \\ c. In order to analyze the data further we 
have traced in Fig. 4 also the universal curve for a clean 
spin-singlet SC with orbital limited upper critical field 
B^^^iO) = 0.72 X Tc \dBc2/dT\T, (WHH model [11]), 
noting that for dirty limit system the prefactor would 
reduce to 0.69. Clearly, the data deviate from the stan- 
dard spin-singlet behavior. Next, we consider the sup- 
pression of the spin-singlet state by paramagnetic limit- 
ing [2^, [2^ . The Pauli limiting field in the case of weak 
couphng is given by 5^(0) = 1.86xrc- For Cuq. 3612.1 Sca 
5^(0) = 6.2 T at ambient pressure. When both orbital 
and spin limiting fields are present, the resulting critical 
field is Bc2(0) = B°^''(0)/Vl% a^^ith the Maki param- 
eter a = V2B°^''(0)/B^(0) [2II23. For B \\ ah (c) we 
calculate a = 1.1 (0.34) and -6^2(0) = 3.3 (1.4) T. These 
values of i?c2(0) are much lower than the experimental 
values (see Fig. 3) and we conclude the effect of Pauli 
limiting is absent. In general, by including the effect 
of paramagnetic limiting the overall critical field is re- 
duced to below the universal spin-singlet values 23l 24|. 
Thus, the fact that our Bc2 data are well above even 
these universal values points to an absence of Pauli lim- 
iting, and is a strong argument in favor of spin-triplet 
SC. The Pauli paramagnetic effect suppresses spin-singlet 
Cooper pairing, as well as the Lz = triplet compo- 
nent (|ti)+ lit))/V2, while the equal-spin pairing (ESP) 
states Itt) and IH) with — I and — —1, respec- 
tively, are stabilized in a high magnetic field. Exemplary 
SCs where Pauli limiting is absent are the spin-triplet SC 
ferromagnets URhCe ^ and UCoGe 

Next we consider the role of anisotropy of the crystal 
structure. Calculations show that for layered SCs, for 
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B parallel to the layers, B°2^ is reduced and the criti- 
cal field can exceed the values of the WHH model [soj . 
Cu2:Bi2Se3 is a layered compound with a moderate 
anisotropy 7"" = 2.9. In this respect it is interesting to 
compare to other layered SCs, like alkali intercalates of 
the semiconductor M0S2 [3l[, which have 7°"-values in 
the range 3.2 — 6.7. A striking experimental property of 
these layered SCs is a pronounced upward curvature of 
Bc2{T) for B parallel to the layers for T < Tc due to 
dimensional cross-over, which is also a salient feature of 
model calculations |30[- However, an upward curvature 
is not observed in Cu2:Bi2Se3. In more detailed theoreti- 
cal work the anisotropy of both the Fermi surface and the 
SC pairing interaction has been incorporated [32l] . Under 
certain conditions this can give rise to deviations above 
the WHH curve as seen in Fig. 4 above. The bulk con- 
duction band Fermi surface of the parent material n-type 
Bi2Se3 is an ellipsoid of revolution along the fee axis with 
trigonal warping [ssj . Assuming Cua;Bi2Se3 has a simi- 
larly shaped Fermi surface, the anisotropy would result 
in a different functional dependence of Bc2{T) for _B || ab 
and i? II c. On the qualitative level this is at variance 
with the universal b*{t) reported in Fig. 4. 

Finally, we compare the Bc2 {T) data with upper criti- 
cal field calculations for a p-wave SC [sj] ■ For an isotropic 
p-wave interaction the polar state (which applies for a lin- 
ear combination of both ESP components) has the high- 
est critical field for all directions of the magnetic field. In 
Fig. 4 we compare the Bc2{T) data with the polar state 
model function. This time the data lie below the model 
curve, but most importantly, the temperature variation 
itself is in agreement with the model, as illustrated by 
the solid black curve in Fig. 4. We have also considered 
a scaled WHH curve, but it fits the data much less well: 
increasing 6*(0) by e.g. 10 % to match the experimental 
value, results in an overall curvature of h*{t) in disaccord 
with the data. In general, topological SC involves all 
three components of the triplet state with a full gap in 
zero field In the case of CuxBi2Se3 model calcula- 

tions indicate triplet pairing is possibly restricted to 
the L2 = component. In an applied magnetic field we 
expect a phase transition or cross-over to a polar state to 
occur. Clearly, more theoretical work is needed on topo- 
logical superconductors in a magnetic field to settle the 
issue of Bc2. 

In summary, we have investigated the pressure varia- 
tion of the superconducting phase induced by Cu interca- 
lation of the topological insulator Bi2Se3. Superconduc- 
tivity is robust and by extrapolating T(,{p) appears to 
vanish at the high critical pressure of Pc = 6.3 GPa. The 
metallic behavior is gradually lost under pressure. The 
upper-critical field i?c2 data under pressure collapse onto 
a single universal curve, which differs from the standard 
curve of a weak-coupling, orbital-limited, spin-singlet su- 
perconductor. The absence of Pauli limiting, the suffi- 
ciently large mean free path, and the polar-state temper- 



ature variation of i?c2 data, point to Cu2;Bi2.iSe3 as a 
p-wave superconductor. 
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